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^ -V ^ Abstract. We fix n and say a square in the two-dimensional grid indexed by 

Qi (x, y) has color c ii x + y = c (mod n). A ribbon tile of order n is a connected 

' polyomino containing exactly one square of each color. We show that the 

^^, ' set of order-n ribbon tilings of a simply connected region R is in one-to-one 

correspondence with a set of height functions from the vertices of R to Z" 
satisfying certain difference restrictions. It is also in one-to-one correspondence 
r'~S ' with the set of acyclic orientations of a certain partially oriented graph. 

' Using these facts, we describe a linear (in the area of R) algorithm for de- 

termining whether R can be tiled with ribbon tiles of order n and producing 
(-H , such a tiling when one exists. We also resolve a conjecture of Pak by showing 

-4-J ' that any pair of order-n ribbon tilings of R can be connected by a sequence 

^ ' of local replacement moves. Some of our results are generalizations of known 

results for order-2 ribbon tilings (a.k.a. domino tilings). We also discuss appli- 
cations of multidimensional height functions to a broader class of polyomino 
tiling problems. 
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1. Introduction 



1.1. Ribbon tilings. A square in the two-dimensional grid indexed by integers 
(a;, y) has color cif x + y = c (mod n). Two squares in the Z^ lattice are adjacent if 
f^ . they share an edge. A ribbon tile of order n is a connected set of squares containing 

"^ I exactly one square of each color. (See Figure ||.) A ribbon tile can also be defined 

'4^ . as a connected sequence of n squares, each of which comes directly above or to the 

2 ] right of its predecessor. A region R is any connected, finite subset of the squares 

of Z^. Two order-n ribbon tilings of R are connected by a local replacement move 
(or simply a local move) if one can be obtained from the other by removing a single 
• rH , pair of ribbon tiles and adding another pair in its place. 

I^% ' Labeled ribbon tilings of Young tableaux (also called rim hook tableaux) have 

^ . connections with symmetric function theory, symmetric group representations, and 

monochromatic increasing subsequences in colored permutations. (See ||l|, 0, g. 
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20 1, and 12^.) Pak was the first to consider ribbon tilings of more general regions 



A ribbon tile of order n has one of 2"~^ possible shapes. (See Figure ^.) If e 
represents one of the 2"^^ shapes of a ribbon tile of size n and a is an order-n 
ribbon tiling of a simply connected region, we denote by a^ (a) the number of tiles 
in a of shape e. 

In [ [L5[ , Pak coined the term tile invariant to describe linear combinations of the 
ae(a) whose values were the same for all tilings a of a region R. Pak showed that if 
it were proved that any two order-n ribbon tilings of a region R could be connected 
by a sequence of local replacement moves, then certain tile invariants of R could be 
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Figure 1. Eight shapes for order-4 ribbon tiles 



deduced as a trivial consequence [|l5| . Pak used a version of the Schensted algorithm 
for rim hook tableaux — a correspondence between rim hook tableau and n-tuples 
of smaller tableaux (known by Nakayama and Robinson ^^ g, rediscovered by 
Stanton and White |2q| ) — to prove that any two order-n ribbon tilings of a Young 
diagram could be connected by a sequence of local moves, and that his tile invariants 
must therefore hold for Young-tableau-shaped regions. Pak used other techniques 
to extend the tile invariants to row-convex regions and conjectured that both local 
move connectedness and (consequently) the tile invariants could be extended to 
arbitrary simply connected regions. 

Without proving local move connectedness, Muchnik and Pak [Q and Moore 
and Pak [Q extended Pak's tile invariant results to (respectively) order-4 and 
order-n ribbon tilings of general simply connected regions. Moore and Pak [nSl 
reiterated Pak's conjecture about local move connectedness for general simply con- 
nected regions and cited private communication from Thurston suggesting that the 
conjecture might be proved if one could find a correspondence between ribbon tilings 
and so-called height functions mapping the vertices of R to Z". Using Thurston's 
suggestion and a connection with acyclic orientations, we prove the conjecture. 

Theorem 1.1. Any two order-n ribbon tilings of a simply connected region R can 
be connected by a sequence of local replacement moves. 

Here, we say R is sim,ply connected ii s ^ R implies that there is no connected 
cycle of squares in R that encircles s with a non-zero winding number; unless 
otherwise specified, we will always assume R is simply connected. 

Our main ribbon tiling result relates the order-n ribbon tilings of a given region 
R to the acyclic orientations extending a certain partial orientation of a graph Gr. 
Pak's original paper |15[ on ribbon tilings used a related construction (based on 
earlier results in S, [ ^0| , and [E5|) to prove local replacement move connectedness 
for ribbon tilings of Young-tableaux. Wc will not use the results of these papers in 
our presentation. 
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Figure 2. A ribbon tiling of order 6, rotated 45 degrees counterclockwise 

Before we can state the main theorem, we need to define several terms. The 
first is a "left of relation, denoted -<, defined for both tiles and squares. In our 
pictures, we generally rotate the set of tiles forty-five degrees so that x + y increases 
as we move upwards. (See Figure ^) Let Sx.y be the square at position {x,y). 

We say Sx.y is directly left of Sx'y ii x + y = x' + y' and x < x' . We say Sx,y is 
indirectly leftoi Sx\y' if \{x+y) — {x' +y')\ — 1 and x < x' and y > y' . (Equivalently, 
Sx^y is indirectly left of Sx' ^y' if it is directly left of one of the four neighbors of Sx\yi ■) 
If Sx,y is directly or indirectly left of Sx',y', we write Sx^y -< Sx'^y' ■ 

If i is a tile and s is a square, we write s ^ t ii s ^ s' for some s' £ t, and i -< s if 
s' ^ s for some s' £ t. If ti and ^2 are two tiles in a, we write ii -< t2 if there exist 
si G ti and S2 € ^2 with si ^ S2. It is easy to verify that if ti and t2 are disjoint 
ribbon tiles, we cannot have both ti -< t2 and f2 ^ ti- 

We will also say a square Sx.y is higher {lower) than Sx'^y' if a; + y > x' + y' 
{x + y < x' + y'). We say x + y is a square of level c if x + y — c and a tile t has 
level c if c is the level of its lowest square. (A square's color is equal to its level 
modulo n.) An edge has type {i,i + I) ii it separates squares of color i and i + 1, 
and a tile has type {i,i + 1) if its lowest and highest squares have color i + 1 and 
i respectively (when defining types, we always add modulo n). Also, a and b are 
comparable by -< if either a ~< b or b ^ a. 

Given a tiling a, we define tc^i{a) to be the ith tile in a, from left to right, of 
level c (when such a tile exists). Clearly tc.i{ci) -< ic.i+i(a) whenever these two tiles 
are defined. 

Let Bn be the boundary of R, that is, the set of squares that are not contained 
in R but are adjacent to at least one square in R. We define the oriented graph 
Gr{ci) to be the graph whose vertices are the tiles in a and the boundary squares 
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of R. Two such vertices a and b are adjacent if they are comparable under ~<, and 
we say this edge is oriented from a to 6 if a -< 6. 

Some aspects of Gala) are independent of a; using height functions, we will 
prove the following result: 

Lemma 1.2. Any two tilings a and f3 of R necessarily contain the same number 
of tiles of each level. Thus, tc.iio) is defined if and only if tc^i{fi) is defined. Fur- 
thermore, if c = d (mod n), then tc,i{a) ~< tdjia) if and only if tc.i{l3) -< td,j{P). 
Also, if b G Br, then tc,i{a) -<b if and only iftc,i{P) -< b . Similarly, b -< tc,i{Q.) if 
and only if b ^ tc,i{P). 

In other words, the map sending tc.iict) into tc,i{f3) is a canonical bijection be- 
tween the tiles of a and those of /3, and this bijection preserves the relation ~< 
between tiles of the same type and between tiles and boundary squares. 

We define a tiling-independent partially oriented graph Gu as follows: all bound- 
ary squares in Br are vertices of Gr and we include a vertex called tc,i in Gr if the 
tile tc,i{a) is defined for some tiling a (and hence all such tilings) of R. An edge is 
contained in Gr if and only if the corresponding edge is included in GR{a) for some 
tiling a (and hence all such tilings). (Thus, ic,i and tdj are adjacent in Gr if and 
only if |c — (i| < n; ic,i and s^^y are adjacent in Gr if and only if c— 1 < x + y < c+n; 
and Sx,y and Sx',y' are adjacent in Gr if and only \{x + y) — (x' + y')\ < 1.) 

Given a tiling a, we can now think of GR{a) as an orientation on the graph Gr 
(i.e., a way of assigning a direction to each edge of Gr) in the obvious way: that is, 
an edge {tc,i,td.j) is oriented from tc,i iotd.j inGR{a) if and only if tc,i (a) -< td,j{a). 
We define the orientation similarly for edges involving two boundary squares or one 
tile and one boundary square. 

An edge involving two tiles of different types (i.e., some ic.i and tdj with c ^ d 
(mod n)) is called a free edge of Gr. All other edges are called forced edges of 



Gr. By Lemma 1.2, a forced edge (a, 5) will be oriented in the same direction as 
GR{a) for every tiling a. Thus, we can think of Gr, as being endowed with a partial 
orientation: its forced edges are all oriented. We say GR{a) extends this partial 
orientation by assigning an orientation to each free edge of Gr. 

Two orientations differ by an edge reversal if they agree on all but exactly one 
edge. The following is our main structure theorem about the set of order-n ribbon 
tilings of R. 

Theorem 1.3. If R admits at least one order-n ribbon tiling, then there is a one- 
to-one correspondence between the order-n ribbon tilings of R and the acyclic ori- 
entations extending the partial orientation of Gr. Furthermore, local replacement 
moves in the space of tilings correspond to edge reversals in the space of orienta- 
tions. To be precise, we make four assertions: 

1. For every order-n ribbon tiling a of R, GR{a) is an acyclic orientation ex- 
tending the partial orientation of Gr. 

2. If a and P are order-n ribbon tilings then a = f3 if and only if GR{a) — Gr{P). 

3. Every acyclic orientation A of Gr that extends the partial orientation of Gr 
is equal to GR{a) for some tiling a. 

4. Two order-n ribbon tilings a and (3 differ by a local replacement move if and 
only ifGR^a) and Gr{P) differ by an edge reversal. 



Given Theorem 1.3, we will now deduce the connectedness of ribbon tilings under 



local replacement moves by citing a corresponding result about acyclic orientations. 
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Define the distance d{A, B) between orientations A and B oia. graph G to be the 
number of edges on which they differ. Then it is well known that A can be trans- 
formed into B with d{A^ B) edge reversals in such a way that each intermediate step 
is also an acyclic orientation. (See g for a more general result.) Clearly, each inter- 
mediate step agrees with A and B on all edges on which A and B agree. It follows 
that if a and (5 are tihngs, then the intermediate steps of a length d{G r^o) , G r{P)) 
path of acyclic orientations connecting Gn^a) and Gr{[3) must all be extensions of 
the partial orientation of Gr. Thus, each of these corresponds to an order-n ribbon 
tiling of R. Assuming Theorem 1.3, we have proved Theorem nTl 



See [pi, p3, p^l and the references therein for more about the structure of 
the space of acyclic orientations. We will also derive, as another consequence of 
Theorem p^.3| , an existence algorithm. 

Theorem 1.4. There is a linear-time (i.e., linear in the number of squares of R) 
algorithm for determining whether there exists an order-n ribbon tiling of R and 
producing such a tiling when one exists. 

We begin in Section g by reviewing the height function theory of Conway- 
Lagarias and Thurston and showing that it leads naturally to a construction of 
abelian height functions for a variety of tiling problems. In Section 0, we apply this 
theory specifically to the ribbon tiling problem and prove our three main results: 



Lemma 1.2, Theorem 1.3, and Theorem 1.4. In SectionM we discuss generalizations 



of our abelian height function constructions to tilings that do not appear to have 
the same acyclic orientation characterization that ribbon tilings have. Finally, in 
Section we present a number of open problems in the theory of ribbon tilings and 
general abelian height functions. 

2. Abelian Conway-Lagarias Thurston Height Functions 

Our height function construction for ribbon tilings makes use of a general tech- 
nique developed by Conway and Lagarias [Q and Thurston [^ for analyzing tiling 
problems. A polyomino is a finite, connected subset of the squares of Z^ whose 
complement is also connected. For simplicity, we will review the theory only for 
polyomino tilings, although similar techniques apply if one replaces the squares of 
the grid with the faces of any planar graph (e.g., the hexagonal lattice). (More 
general expositions — which invoke the language of Cayley complex homology — can 
be found in [|l), |2|, and §.) 

A tile set T is a (usually infinite) set of finite, simply-connected subsets of the 
squares of Z^. Usually, we assume T is translation invariant (i.e., t G T implies that 
{sx,y\Sx+xo.y+yo ^ t} G T for all xq, j/o 6 Z) but this need not be the case. The first 
step in the Conway-Lagarias construction is to choose a map </> from the oriented 
edges of the plane into a group G in such a way that (j){a) = —(j){a~^) for all oriented 
edges a; here a~^ is the same edge as a with opposite orientation, and (since we will 
eventually restrict our attention to abelian groups) we write the group operation 
additively. We also require that the sum of (j){e) over the clockwise-oriented edges e 
on the boundary of any tile in T (added in clockwise order) be equal to the identity 
in G. These requirements are called tile relations. Once the group is chosen, given 
any tiling a, we can construct a height function ha mapping at least some of the 
vertices of R into G as follows. 

First, choose a reference vertex vq on the boundary of R and require that ha{vo) 
be the identity. Then if v is any other vertex of R that is on the boundary of a 
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tile, ha{v) is the ordered product of the oriented edges along a path from vq to v 
that does not cross any tiles. (If the tiles are simply connected, such a path will 
always exist.) Using the tile relations, it is not hard to prove that the value ha{v) 
is independent of the path chosen. 

We can also always assume that the group G is generated by the values of 
(j){e) where e is an edge in the Z^ lattice. (Otherwise, replace G by the subgroup 
generated by these values.) One can obtain the largest possible group G that is 
generated by (j){e) as follows: Let G be the group generated by the oriented edges of 
7? modulo the relations ^ e = for clockwise sums of clockwise-oriented boundary 
edges of a tile and e + e^^ == for all e. Then let ^(e) be the image of the edge e 
in the group G. 

The height function constructed in this way is in a sense "maximally informative" 
since any other height function is a quotient of this one. Unfortunately, this group 
is unwieldy in practice and can be difficult to analyze. On the other hand, if one 
merely seeks the maximally informative abelian group, the analysis becomes much 
simpler. 

Suppose G is an additive abelian group and is a map from the oriented edges 
of Z^ to G. If t is any polyomino tile in 7?, we denote by (/)(t) the sum of the 
clockwise-oriented edges of t. Observe that 

0(i)-5]0(s), 

where the sum is over squares s that are contained in the tile t. 

Now, suppose that p = (wq, wi, . . . , f^) is a non-self-intersecting path of vertices 
in R connecting the reference vertex vq to a vertex Vk on the boundary of some tile 
in a. We say that p crosses a tile f in a if for some < i < fc — 1, the edge {vi, Wi+i) 
lies between two squares in t. If s is a square in R and < i < j < k, we say p 
passes s on the left between vertices i and j if: 

1. The vertices Vi and Vj are on the boundary of a tile t and all edges between 
Vi and Vj on p are interior edges of the tile t. 

2. The square s is contained in the polygon formed by the edges in p between 
Vi and Vj together with the counter-clockwise half of the boundary of t con- 
necting Vi to Vj. 

Note that a long path could pass a square s on the left several times if s belongs 
to a tile that the path crosses several times. 

We define 4'{p) = X]i=o 0(^ii^i+i) ^'^d observe the following: 

Proposition 2.1. The height function value at Vk can be written as 

K{vk) = 0(p) + ^(/'(s), 

where the latter sum is over all squares (counted with multiplicity) that p passes on 
the left. 

To see this, let p' be a path from wq to Vk that is the same as p except that each 
interval that crosses a tile t is replaced by the path around the boundary of t on 
the left side of p. The reader may easily check that (j){p') — 4>{p) is equal to the sum 
of (/)(s) over all squares (counted with multiplicity) that p passes on the left. By 
definition, ha(vk) = 4'{p'), so the proposition follows. 

Since (j)(j)) is independent of the particular tiling a, we have a simple interpre- 
tation of the height function ha{v): it tells us the sum of (j){s) over the squares s 
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Figure 3. This path of vertices (beginning with the lowest bound- 
ary vertex) passes each shaded square on the left once. 



passed on the left by a path from vq to w. For our height function to be maximally 
informative, G must be (or at least contain) the abelian group generated by (j){s) 
subject to only the tiling relations. We refer to this group as the canonical abelian 
height space (also known as the color group) of a set of tiles T. 

We say two squares si and S2 have the same color if it can be deduced from the 
tiling relations that 0(si) = 0(s2)- (In other words, si — S2 is zero in the group 
generated by the squares of Z^ modulo the relations '^g^^. s = for all t G T.) 

This notion is best illustrated using our primary example: ribbon tiling. Let 
Sx.y denote the square at position {x, y) in the lattice, and consider two overlapping 
horizontal ribbon tiles of length n, one of which has leftmost square s^.y and the 
other Sx+i,y. The corresponding tile relations tell us that X]i=o 4>{sx+i,y) ~ and 
m=i 4>{sx+i,y) = 0. And this implies that 4>{sx,y) = 4>{sx+n,y)- Thus, Sx,y and 
Sx+n,y have the same color. 

Repeating similar arguments with other pairs of ribbon tiles that overlap on all 
but one square, one can deduce that squares at {xi,yi) and (a;2,y2) necessarily 
have the same color whenever xi + j/i and X2 + 2/2 are equivalent modulo n. Denote 
by Ci the value of 4'{s,j..y) when a; + y = z (mod n). Then the relation X]r=o ^i ~ ^ 
is enough to guarantee that (j){t) = for any ribbon tile t of length n. Thus, if 
G is the abelian group generated by cq, . . . , e„_i subject to this one relation, G is 
maximally informative. 

At this point, we know the value of (j){s) G G for every square s. However, we 
still have to choose a function (p on the edges of the 1? lattice in such a way that 
(/)(s) has the desired values. This is essentially a linear algebra problem, and it is 
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not hard to see that it always has at least one solution. (We return to this issue in 
the next section.) 

This "maximally informative" abelian height function construction is unique up 
to two obvious transformations. First, change of basis: one can replace G with 
any group G" and the e^ with any elements of G" satisfying the necessary relations. 
Second, translation by a tiling-independent function: one can add to 4'{e) any 
function <f>' satisfying 4>'{s) — for all s. This amounts to adding a function 
(independently of a) to each height function ha- 

In particular applications, we try to choose the basis for G and the values of cf) 
on the edges in 1? in a way that makes the tiling space easy to understand. In 
particular, if T is translation invariant, it is usually convenient for (p to have some 
translational symmetries as well. 



3. Ribbon Tilings 

3.1. Height Functions for Ribbon Tilings. Although we determined in the last 
section that 4){s) — ei when s is a square of color i, we have considerable freedom 
in choosing the values of (f) on the edges in 7? . We will use this freedom to force (/) 
to have certain symmetries. 

First, we would like (j) to be invariant under color- preserving translations of 1?] 
this amounts to the requirement that 4>{e) be the same whenever e is a vertical 
(or similarly, horizontal) edge of type (i,i -f 1) oriented with color z -I- 1 on the 
left. We would also like ((> at such an edge to have the same value — call it e^.i+i — 
independently of whether the edge is horizontal or vertical. 

Let eo,i, ei^2, • • ■ , e„„i.o be basis elements of a height space G = Z". Summing 
around a square of color i tells us that e^ must be equal to 2ei^i+i — 2ei_i.i. Since 
this clearly implies X]r=o ^« ~ ^^ ^^ need not add any additional relations to G. 
(Although the canonical abelian height function space is the subgroup of G gener- 
ated by the values e^ — 2ei_i+i — 2ei_i.i, it will be notationally convenient to define 
(j) and our height functions in this larger group.) We can now describe our height 
function ha for a tiling a with the following rules: 

1. Let /iq(wo) = for all tilings a of R, where uo is a reference vertex on the 
boundary of R. 

2. When moving from one vertex to a neighbor along an edge separating i and 
i -\- 1 color squares with color i + 1 on the left, the height changes by e^^i+i if 
no tile of a is crossed. 

3. When moving from one vertex to a neighbor along an edge separating i and 
i + 1 color squares with color i -I- 1 on the left, the height changes by 2ejj+i — 
Ci^i+i if the edge crosses a tile in a of type (j, j + 1). 

We write /i^*+^ to mean the e^.i+i component of the height function ha- Of 
course, there is some redundancy in using an n dimensional height space instead 
of the canonical n — 1 dimensional height space; as currently defined, Y11^=q ^q*^^ 
is a function that does not depend on the tiling a. In particular, if n = 2, the 
projection sending aeo,i -I- hei^ to the integer a — h transforms our two-dimensional 
ha into the standard one-dimensional domino-tiling height function without any 
loss of information. 
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3.2. Interpretations of Ribbon Tiling Height Functions. We can understand 
what the height function means by looking at its behavior along diagonals. We say 
p is an (i, J + 1) diagonal if: 

1. The path p — vq, . . . , uj. is a connected left-to-right path made solely of edges 
of type (i, i + 1) that are incident to squares in R. 

2. Both vq and Vk are on the boundary of R and the extending edges (u_i,wo) 
and (ufc, Vk+i) of type (i, z -t- 1) are not incident to squares in R. 

Note that as one moves from left to right along an {i.i + 1) diagonal, ft-^'-'^^ 
increases by two if we cross a a tile of type {j,j -\- 1) and zero otherwise. Thus, in 
the diagonal described above, h^^^~^^{vi) is equal to h^a^~^^{vo) plus twice the total 
number of times the path vq, ■ ■ ■ ,Vi crosses a tile of type {j,j + 1). In particular, 
since the value of h^ on the boundary of R is independent of a, we have the 
following. 

Lemma 3.1. If j ^ i, then the total number of tiles of type (j, j + 1) crossed by the 
portion of an (i,i + l) diagonal between two consecutive boundary vertices vi and V2 
is equal to ■^{h-'^-^^^{v2) — h^g^^~^^{vi)). In particular, since Vi and V2 are boundary 
vertices, this number is independent of a. 

How docs h"^^ change along an {i, i + 1) diagonal? We know that moving from 
left to right along the diagonal, it decreases by one when we cross a tile, increases 
by one otherwise. Also, inspection shows that there can never be tiles crossing both 
of two adjacent edges on a diagonal. Thus, as one moves along the diagonal from 
left to right, the value of /i^*"''^ alternates between sequences of decreasing steps (of 
length at most one) and sequences of increasing steps (of length one or more) . We 
say Va is a record vertex if a > 1, and /i^*~''^(ub) < h'^a'^^iva) whenever < 6 < a. 
For notational convenience, we will also consider Vk+i (the vertex just right of the 
vertex Vk along an edge of type (i, z -I- 1)) to be a record vertex. (See Figure ^) 

Since each decreasing sequence in /i^*"*"^ along p has length at most one, it is 
not hard to check that if 2 < a < fc, then Va is a record vertex if and only if 
h'^^'^'^iva) > ^ji*"'""^(fa--i) > hi;'-+^{va-2)- Notc that the first edge (t'o,wi) and last 
edge (ffe-i, Vk) of a diagonal are always boundary edges of R, and hence represent 
increasing steps. In particular, vi is always a record vertex. 

What does it mean for Va to be a record vertex? First of all, it means that the 
two adjacent edges {va-2,Va-i) and {va-i,Va) are both not crossed by tiles; this in 
turn implies that the square s incident to both of them (which we will call a record 
square) is either a square outside of i? or a square in R belonging to a tile of type 
(i,z + l). Since each of the integers /i^'+^('(;o)-|-l, /i^'"'"H^o) + 2, . . . , /i^*"'""^(ufe) must 
occur as a record value exactly once, the number of record vertices on the diagonal 
(not counting Vk+i) is equal to /i^*"'""^(wfe) — h^a~^^{vo), and is hence independent of 
a. 

Now, if s is a record square outside of R and Va is the corresponding record 
vertex, the two diagonal edges preceding Va must both be boundary edges of R 
(unless one is an edge outside of i?); hence, s will be a record square for every 
tiling of i?. If a 7^ fc + 1, then the value ha{va) is also independent of the tiling. 
Whenever s is a record square outside of R, the corresponding record vertex is 
called a boundary record vertex; in particular, vi and Vk+i are boundary record 
vertices. 

Between any two consecutive boundary record vertices Va and vt, there are ex- 
actly hl;^~^^{vb) — /iQ*'^^(wa) — 1 record vertices whose record squares are squares 
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Figure 4. Let p be the diagonal of type (2, 0) that separates the 
shaded squares from the unshaded squares. If vq is the leftmost 
corner of the figure, then h^^{vo) — 0. As one moves from left 
to right along the diagonal, ft.^'° assumes, in order, the values 
0, 1, 2, 3, 4, 5, 4, 5, 4, 5, 6, 7, 8, 7, 8, 7, 8, 7, 8, 9, 10, 9, 10. The dots in- 
dicate record vertices — one for the first occurrence of each number 
between 1 and 10 and one for the first vertex outside of R\ the 
square immediately left of each record vertex is either a boundary 
square or a square in a tile of type (2,0). The relation -< induces 
a total ordering on this set of four boundary squares and seven 
type (2,0) tiles. Every pair of elements in this set corresponds to 
a forced edge of G_R(a). 



of tiles of type (i, i + 1). Whether the tile lies above or below the diagonal clearly 
depends on the parity of a; since /i*'*+^ alternates parity along the diagonal, this 
in turn depends on the parity of the /i^*+^ at the record vertex Va- We call v a 
lower record vertex if the corresponding record square lies below the diagonal and 
an upper record vertex if it lies above the diagonal. Since for each a, the record 
values h'^a^^iyo) + ^■•^'^a'^^i'^o) + 2,... ^h^a'^^i'^k) must be attained in the same 
order, we have the following lemma: 

Lemma 3.2. Suppose q is the portion of a diagonal p of type {i,i + 1) between 
two consecutive boundary record vertices Va and Vh. Then the number tiles of type 
{i,i + 1) with squares adjacent to q that lie above or below the diagonal is equal 
to hl^^~^^{vb) — /i^*^^(wa) — 1. In particular, this value is the same for every tiling 
of R. Furthermore, the order in which the entire diagonal p encounters boundary 
record vertices, lower record vertices, and upper record vertices is also the same for 
every tiling of R and can be determined from the value of the height function on 
the boundary vertices in p. 

Denote by Sa{p) the set of squares adjacent to p that are either members of tiles 
of type {i,i + 1) or squares in Br. (See Figure ^.) Note that if Si and S2 are two 
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record squares then si ^ S2 if and only if the record vertex for si occurs before 
that of 32 in the diagonaL Thus, we can restate the previous lemma: 

Lemma 3.3. Given any two tilings a and a' , Sa(j>) and Sa'{p) contain the same 
number of boundary squares, squares of color i in type {i,i + 1) tiles, and squares 
of color i + 1 in type {i,i + 1) tiles. The relation -< induces a total ordering on 
the squares in Sa{p) as well as those in Sa'{p); if the mth square in the ordering 
of Sa{p) is a boundary square (or a square of color i in a type (i,i + 1) tile, or a 
square of color i + 1 in a type {i,i + 1) tile), the same is true of Sa> (p). 

Now we are ready to prove Lemma n^, which we restate here for convenience. 

Lemma 3.4. Any two tilings a and f3 of R necessarily contain the same number 
of tiles of each level. Thus, tc,i{a) is defined if and only iftc^i{f3) is defined. Fur- 
thermore, if c = d (mod n), then tc.i{a) ~< td.j{a) if and only iftc.iiP) -< td.jiP). 
Also, if b (z Br, then tdict) < b if and only if tc,i{l3) -< b . Similarly, b -< ic.i(«) */ 
and only if b ~< tc,i{l3). 

Proof. Let ic,i(a) ^ ^c,2(ck) -<...< tc^kicn) be the tiles of type (i, z + 1) in a whose 
lowest squares lie just above a given diagonal p. It is obvious from the definition 
of tc^i{a) that the ordering of these tiles with respect to each other is independent 
of a. Similarly, let tc-n,i{ct) -< tc~n,2{ct) < . . . ^ tc^n,m{a) be the tiles of type 
{i, i + 1) that lie just below the diagonal p. 

Now, let ai ~< 02 < . . . ~< Or be the totally ordered set of a^Zthe boundary squares 
and all of the k + m tiles that are incident to p. Since we know that the ordering 
restricted to the set of tiles of form tc,i{a) (or the set of tiles of form tc-nA^), or 



the set of boundary squares) is independent of a, it follows from Lemma 3.3 that 
the total ordering ai ~< 02 ~< . . . < Or is independent of a. 

Since any two tiles of the same type or any tile and any square that are com- 
parable by -< must be incident to some common diagonal p, Lemma ^.3| covers all 
possible cases in which the -< relation occurs between tiles of the same type or 
between tiles and boundary squares. D 

The reader may recall from the introduction that we used this lemma used to 
define the partially oriented graph Gn and to justify our interpretation of Gru as 
an extension of the partial orientation of Gr. This will be the topic of the next 
section. 

3.3. The Acyclic Orientation Correspondence. The purpose of this section 



is to prove Theorem 1.3; we will state and prove each of its four parts as a separate 
lemma. 

Lemma 3.5. For every order-n ribbon tiling a of R, Gii{a) is an acyclic orienta- 
tion extending the partial orientation of Gr. 

Proof. We must show that, given a tiling a, there exists no sequence 

a^ < ai ^ 02 < ■ . . < Ok < Oq, 

where each a^ is either a tile or a square in Br. 

Since there can clearly be no cycle containing one or two elements, we may 
assume fc > 2. Now, suppose a diagonal incident to both ag and a^ is also incident 
to some tti with Q < i < k. li Ok < flj, then Oi -< a^+i -< . . . ^ Ok ^ Oi is a shorter 
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cycle. If Qi -< ttk, then a^ -< ao (by transitivity, since Ui, Uk, and ao are on the same 
diagonal and hence totally ordered), and oq ^ ai -<...-< a^ -< ag is a shorter cycle. 
Thus, if ao ^ ai -< a2 -< . . . ^ afc ^ ao is a cycle of minimal length in G/{(a), 
then any infinite diagonal that is incident to both ag and ak can be incident to 
no other tile or square. Denote by Pz the infinite diagonal between square levels 
z and z + 1 and let z' and z" be the lowest and highest values of z for which p^ 
is incident to both ao and Uk- AH tiles and squares other than ao and a^ must lie 
entirely below Pz' or entirely above Pz" ■ Since no tile or square below pz' can be 
comparable with any tile or square above Pz", the sequence ai ^ a2 ^ ■ ■ ■ ^ o-k-i 
must lie entirely above Pz" or entirely below Pz' ■ Suppose without loss of generality 
that the former is the case and that the lowest square in ao is at least as low as the 
lowest square in ak- Then none of the elements a^ below Pz' is comparable with a^, 
and this is a contradiction. D 

Lemma 3.6. // a and (3 are order-n ribbon tilings then a = (3 if and only if 
GR{a) = GR{l3). 

Proof. Obviously, the first statement implies the second. Now, suppose Gii{a) — 
Gn{P). Since these orientations are acyclic, there must exist at least one leftmost 
tile, i.e., there exists a ic,i for which there exists no i^j with tdj{a) -< tc^i{a) or 
td,j{f3) -< tc,i{f3). It follows that tc.i(a) and tc.iiP) must each contain the leftmost 
squares in R in the rows of level c, c+l,...,c + ri — 1; thus, their positions are 
determined and they are equal. Similarly, we can now choose a i^j that is leftmost 
among the remaining set of tiles, and its position is also completely determined. 
Repeating this argument, we conclude that tc,i{a) = tc^i{P) for all c and i for which 
these values are defined. D 

Lemma 3.7. Every acyclic orientation A of Gr that extends the partial orientation 
of Gr can be written as GR{a) for some tiling a. 

In the introduction, we defined Gr by assuming that R admitted at least one 
tiling a and restricting the orientation on GR{a) to the forced edges. However, 
in our later existence proofs, it will be necessary to have a definition of Gr and a 



version of Lemma 3.7 that do not invoke a priori knowledge of the existence of a. 
To do this, we will need some definitions. 

We say a tile vertex ic,i crosses level difc<d<c + n — 1. (Note that if a is 
a tiling, then tc,i crosses level d if and only if ic,j(a) contains a square of level d.) 
We define a partially oriented graph G to be a ribbon tile graph consistent with R 
if the following hold: 

1. The vertices of G come from the set of tile vertices ic,i (with c,i E Z and 
z > 0) and square vertices Sx,y (with x,y gZ). 

2. Two tile vertices tc,i and t^j are adjacent if and only if |c — dj < n. A tile t^^i 
and square s^^y are adjacent if and only iic— l<x + y<c + n. Two squares 
Sx,y and Sx',y> are adjacent if and only if |(a; + y) ~ (x' + y')\ < 1. 

3. The graph G is endowed with an orientation on its forced edges (i.e. all 
edges involving a square or involving two tile vertices tc,i and tdj with c = d 
(mod n)). No other edges of G are oriented. 

4. If si and S2 are square vertices, then the edge (si, S2) is oriented toward S2 if 
and only if si -< S2- Also, whenever i > and ic,i+i is a tile vertex, tc,i is also 
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a tile vertex. An edge of the form {tc^i,tcj) is oriented toward tcj whenever 
i < j. 

5. The set of square vertices of G is precisely the set of boundary squares of R. 

6. The total number of tiles crossing a given level c is precisely the number of 
squares in R of level c. 

7. For every integer c and every pair of boundary squares 61 and 62 of levels 
c — 1 , c, or c + 1 satisfying 61 ^ 62 , the number of tile vertices tdj in Gh that 
cross level c and satisfy 61 -< t^j -< 62 is precisely the number of squares s in 
R of level c satisfying 61 -< s ^ 62 . 

Only the last three conditions depend on R; we will refer to these as the con- 
sistency conditions. Although it involves overloading the symbol -<, we will write 
a < b whenever the edge (a, b) is a forced edge of G oriented toward h. If A is 
an orientation extending the partial orientation on the forced edges of G, we write 
a -<A b if the edge (o, b) in G is oriented from a to 6 in A. 

Note that up to this point, we have defined -< only on actual tiles, not on the 
tile vertices tc^i, which are abstract place holders for tiles. However, our expanded 
use of -< is consistent in the following sense: if Gij(a) is the ribbon tiling graph 
of a tiling a and the edge {tc.i,td,j) is forced, then tc,i{a) -< tdjiP) if and only if 
tc,i -< td,j in Gr. Similarly, if ^ = Gfl(a) and the edge is (ic,i, ^dj) is not necessarily 
forced, then ic,i(a) ^ id.j{ct) if and only if tc,i <a td,j- Similar results apply to 
edges involving squares. 

If R admits a tiling a then the graph Gr as defined in the introduction clearly 
satisfies the above list of requirements and is consistent with R. (In fact, Gr is 
the only ribbon tile graph consistent wi th R in this case.) Thus, the following is a 



slightly more general version of Lemma 3.7 



Lemma 3.8. Suppose Gr is any ribbon tile graph that is consistent with R. Then 
every acyclic orientation A of Gr that extends the partial orientation on the forced 
edges of Gr is equal to GR^a) for some tiling a. (In particular, if R does not admit 
an order-n ribbon tiling, then no such Gr and A exist.) 

Proof. Suppose t — tc.i is a, leftmost tile vertex in Gr. As in the previous proof, 
we will choose the leftmost squares of each of the levels c, c+1,... ,c + n — l and 
declare them to be our first tile tc,i(a). That there exists at least one square in 
R on each of these levels follows from the consistency conditions of Gr. However, 
we must verify that this set of squares is actually a ribbon tile. Since it contains 
one square on each of levels c, c + 1 , . . . , c + n — 1 , it is sufficient to show that it is 
connected. 

Suppose otherwise; then we can assume without loss of generality that there 
exists some i such that the leftmost square Si in the row of level c + i is left of the 
square s^+i in the row of level c + i + 1 but that these two squares are not adjacent. 
Now, let Sr (s/) be the square of level c + i + 1 that is right of (left of) and adjacent 
to Si. Both Sr and Si must be squares in Br, since they are adjacent to Si and they 
are left of s^+i, the leftmost square of their level in R. 

Now, by the consistency conditions, there is exactly one tile t' containing a square 
on level c + i such that si ^ t' ^ s,- ; every other tile with a square on level c + i 
must be right of s^. If t ^ i', then t' < Sr < t ^ t' , violating acyclicity. Hence 
t = t' . However, since Sr is left of every square of level c + i + 1, one may deduce 
from the consistency conditions that every tile vertex t of Gr that crosses level 
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c + i + I must satisfy Sr ^ t. Thus t ^ Sr ^ t, a, contradiction, and it follows that 
the squares chosen do in fact form a tile. 

Now we form a new region R' by removing the squares in this tile from R. We 
form a new graph Gr/ by removing fc,i from the vertex set of Gn — along with 
all boundary squares of R that are not boundary squares of R' — and adding all 
the boundary squares of R' that were not boundary squares of R. (These are 
necessarily members of our removed tile.) We define the edges of Gr' and their 
orientations in a new orientation A' in the most natural possible way. Two vertices 
of G/}/ are adjacent whenever the corresponding edges in Gji were adjacent, and A' 
is the orientation induced by A on such edges. Furthermore, if s is one of the new 
boundary squares and is in a row c, we declare it to be adjacent in Gri to every 
tile t that crosses rows c — 1, c, or c + 1 (with s ^ t) and every boundary square 
So of those levels (with s ^ sq if and only if s is directly or indirectly left of sq and 
sq ^ s when the opposite is true). Finally, if tcj was a vertex of Gr with j > 1, 
we relabel it tcj-i in Grl 

It follows easily from the definition that the Gri thus defined — together with the 
partial orientation on forced edges induced by A' — is a ribbon tile graph consistent 



with R'. If we can further prove that A' is necessarily acyclic, then Lemma 3.8 will 
follow by induction. 

Suppose A' has a cycle; then this cycle must include a square s of Br/ that was a 
member of the tile t — otherwise the cycle would also be a cycle in A. Also, although 
the cycle may contain multiple squares from the removed tile, it must contain at 
least one element that is not a square of t, since the squares do not themselves 
contain any cycles. Suppose the cycle contains a sequence si ^ . . . ^ s^ of former 
squares of t but that the elements a immediately preceding this sequence and b 
immediately following it are not squares of t. Then we must have had a ^ t ^ 6 in 
A, and A too must have had a cycle. D 

Lemma 3.9. Two order-n ribbon tilings a and [3 differ by a local replacement move 
if and only if GR{a) and Gr{P) differ by an edge reversal. 

Proof. First of all, we can determine whether a local replacement move involving 
tiles tc,i{o:) and tdj{a) can occur by applying the theory we have already developed 
to R' , the region tiled by just these two tiles. Now, if R' admits two distinct tilings, 
then Gri must admit two distinct orientations that extend the orientation on the 
forced edges. In particular, Gr' must contain at least one free edge, and that can 
only be the edge (ic,i, td,j)- For this edge to be free, we must have < |c — dj < ti. 
Without loss of generality, suppose c < d. Notice that if(i<e<c + 7i— 1, then 
R' contains two squares of level e. Inspection shows that if tile ic,i(a) contains the 
leftmost square of one of these same-level pairs, it must contain the leftmost square 
of every such pair. (See Figure o.) Moreover, if a tile contains one square on a 
row higher than c + n — 1 (lower than d) , it must contain every square with this 
property. Thus, a local replacement move must consist of swapping the squares in 
each same-level pair. Such a swap produces a new pair of ribbon tiles if and only 
if the single square on row d— 1 (resp., c + n) is adjacent to both of the squares in 
row d (resp., c + n — 1). Inspection shows that whenever there is no square s with 
either tc^i{a) < s < tdj{a) or tdj{a) -< s ^ tc^i{a), then a local move can in fact 
occur. 
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Figure 5. Two examples of local replacement moves 

Now, if f3 is the tiling obtained by applying such a move to a, and tc.i{a) -< t^j (a) 
then inspection shows that we must have tdj{P) -< tc,i{f3); it is also easy check that 
if a is any tile or square outside of R', then a -< tc,i{cx) if and only if a -< tc^i{P)- 
Similarly, tc,iia) ^ a if and only if tc^i{f3) -< a. It follows that if a and (3 differ 
by a local replacement move, Guia) and Gb{I3) do in fact differ by the reversal of 
exactly one edge: the edge {tc,i,tdj)- 

For the converse, note that if the orientation of an edge {tc.i,td,j) of the ori- 
entation Gji{a) can be reversed without creating a cycle, we know there must be 
no vertex a of GR{a) with either tc,j(a) < a < td.j{a) or tdjia) -< a -< tc,i{oi); 
we have already shown that this is a sufficient condition for there to exist a local 
move involving the tiles tc,i{a) and tdj{a), and so Gb.{I3) must be the orientation 
obtained by reversing the edge. D 

3.4. The Existence Algorithm. The main result of this section will be a linear- 
time algorithm for determining whether a simply connected region R has an order-n 
ribbon tiling and constructing such a tiling when one exists. To determine existence, 
we will make use of the well-known fact (which we later prove) that there exists 
an acyclic orientation A of a graph G that extends a given partial orientation of 
G if and only if the par tial orientation is itself acyclic. This implies the following 
corollary of Lemma 3 



Corollary 3.10. Given R, there exists an acyclic partially oriented ribbon tiling 
graph Gr that is consistent with R if and only if R admits an order-n ribbon tiling. 

Using this fact, we divide our algorithm into three steps. 

1. We either produce an acyclic ribbon tile graph Gr that is consistent with R 
or prove that no such graph exists (and hence R admits no order-n ribbon 
tiling) . 

2. If we succeed in producing Gr in the previous step, we describe a total or- 
dering on the vertices of Gr that induces the partial orientation of Gr on the 
forced edges of Gr. 
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3. Given this totally ordering on the vertices of Gr, we describe a tiling of R. 
Now, we know that if there is an acyclic ribbon tile graph Gn that is consistent 



with R, we can use it to generate Gii{a) using the construction of Lemma 3^, and 
the partial orientation on Gr will thus be the restriction of G'ij(a) to the forced 
edges in Gr. Thus, in the construction that follows, we can make use of lemmas 
proved earlier in the paper under the assumption that R had a tiling. If R does 
have a tiling, our assumption will have been correct. If R does not have a tiling, 
these assumptions will lead to contradictions, and these will enable us to conclude 
that no tiling exists. 

Now we begin to construct Gr. Let S{c) be the number of squares of level 
c and let T{c) be the number of tiles whose lowest square is c. Clearly T{c) — 
S{c) — Y^Zi ^(c ~ i), and is thus easily computed for each c in linear time. (If it is 
not the case that T{c) > whenever S{c) ^ and T{c) = whenever S{c) = 0, we 
can conclude that no tiling exists.) For each c for which T{c) ^ 0, we include tc,i 
for 1 < i < T{c) as tile vertices of Gr. We also include all the boundary squares 
of R as square vertices of Gr. This Gr has an edge whenever the corresponding 
tiles/squares are on levels that would make them comparable by <, and we require 
that tc.i ^ tc,i+i whenever both of these vertices are in Gr. 

We would now like to determine an orientation on the forced edges of Gr that 
will make Gr a consistent ribbon tiling graph. We will not store in memory all of 
the edges of Gr (as the number of such edges could be quadratic in the size of R) ; 
instead, we will observe that the partial orientation on Gr induces a total ordering 
ao -< ai ^ . . . ^ ttk on the set of squares and tiles of a given type that are all 
incident to a common diagonal. Since every forced edge can be derived from such 
an ordering on some diagonal, it will thus be sufficient for us to store in memory the 
actual orderings corresponding to each such type/diagonal pair, and this requires 
only a linear amount of memory. 

Given an infinite diagonal p between square rows of levels c — 1 and c and some 
j G {1,2,... ,n — 1}, let Sj^p be the set of tiles in Gr of the form tc-j^i together 
with squares in Gr, on levels c — 1 and c. Since the elements in Sj^p share a common 
diagonal and all tiles in this set are of the same type, each pair of elements in Sj^p is 
connected by a forced edge and hence Gr induces a total ordering on this set. Now, 
we can compute the height function h on the boundary of R, and using h, we can 
deduce this ordering explicitly from Lemma [3. l| as follows. Let k — c—j—1 (mod n). 
Given consecutive boundary vertices vi and V2 along the diagonal, we know that 
the number of tiles of the form ic-i,i (and hence of type {k,k + 1)) that cross p 
between consecutive boundary vertices wi andw2 must be ^{h'''''^^{v2)—h'''''^'^{vi)). 
Clearly, if 5i is a square incident to the diagonal p immediately before vi and 62 is 
a square incident to p immediately after V2, the number of tiles t of the form ic-j,i 
satisfying 61 -< i ^ 62 must be exactly ^(/i'''''+^(w2) - /i'''''+^(wi)). Thus, we can 
compute the number of tiles of Sj^p in between each pair of boundary squares, and 
this completely determines the ordering. 

Similarly, given an infinite diagonal p between square rows c — 1 and c, let iSp_o 
be the set of tiles in Gr of the form tc,i or tc-n,i together with the squares in Gr 
on levels c — 1 and c. Each pair of elements in iSp_o is connected by a forced edge, 
and thus the partial orientation on Gr must induce a total ordering on this set. 



By Lemma B.2 and Lemma 3.3, we can deduce this total ordering from the value 



of the height function on the boundary vertices of R in p. 



RIBBON TILINGS AND MULTIDIMENSIONAL HEIGHT FUNCTIONS 17 




Figure 6. A graphical representation of Gr when n = 3 and R is 
the 9x9 grid. Each of the dots represents one of the 27 tile vertices 
of Gr (square vertices are not shown), and the 15 horizontal lines 
correspond to the 15 levels on which tiles occur. All essential forced 
edges between tiles of the same type are shown with orientations. 
Two vertices in this figure are adjacent in Gr if and only if the 
difference between their levels is at most three. 

Since every forced edge involves two elements along the same diagonal, we can 
compute the orientation of each forced edge by repeating this process for each 
diagonal. Now we check that if Gr is acyclic, it is necessarily a consistent ribbon 
tiling graph. It is enough to check the consistency conditions. By construction, 
Gr contains all boundary squares, and by the definition of T{c), the total number 
of tiles crossing each level c is the number of squares in R of level c. For the final 
consistency condition, suppose c is given and bi and 62 are boundary squares of 
types c — 1, c, or c + 1 with bi -< &2- 

Then let pi be the diagonal above the row of squares of level c and let p2 the 
diagonal below these squares. One easily checks that both pi and p2 contain edges 
in each of &i and 62- It suffices to check the condition under the assumption that 
there exists no boundary square 63 incident to p with 61 -< 63 -( 62, since otherwise 
we can check the condition separately for the pair (61, 63) and the pair (63, 62). Say 
a vertex has level c if it is a corner of squares of type c — 1, c, and c + 1, and let 
q be the portion of the diagonal connecting the last level c vertex of 61 to the first 
level c vertex of 62 ■ 



Then it follows by applying Lemma 3.1 to the diagonal pi (or to p2 if pi has type 



(j, j + 1)) that the number of tile vertices t of type (j, j + 1) crossing c and satisfying 
61 -( t -< 62 must be equal to ^{h^^^+^{v2) - h^^^+\vi)). Write h = ^^'ig ft^'^+i, 
and note that the total number of such tile vertices must then be ^{h{v2) — h{vi)). 
However, this function is tiling-independent and in fact always increases by one as 
we move between squares of color j and j + 1 with j + 1 on the left. It follows that 
if the vertices of the grid are labeled with coordinates (u,w), then h{u,v) will be 
equal to u — v plus a constant. In particular, this implies that ■k{h{v2) — h{vi)) is 
equal to the number of squares on level c directly between vertices vi and V2 , and 
this is easily seen to be the number of level c squares s in i? satisfying 61 -< s -< &2 ■ 
The only remaining step is to determine whether our partial orientation of Gr 
is indeed acyclic. There are many of ways to do this in time linear in the number 
of edges of Gr. We say a forced edge (a, 6) of Gr is essential if a -< & and there 
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exist no forced edges (a, c) and (c, 6) with a ~< c and c -< b — one easily checks 
that the number of essential forced edges of Gr is linear in the size of R and that 
any acyclic orientation of Gn that extends the partial orientation on the essential 
forced edges must extend the partial orientation on all forced edges. (See Figure |[) 
The following algorithm for determining acyclicity (know as the topological sort) is 
well-known, and we will apply it to the graph formed by the vertices of Gr and the 
essential forced edges of Gr. 

Choose a starting vertex oq from the set of vertices of Gr and attempt to con- 
struct an increasing path oq -< ai ^ 02 ■ ■ ■ until we either reach a point we have 
already visited (and we declare that there exists a cycle) or we reach a vertex a^ 
that is not left of any of its neighbors. Clearly, Ok cannot be part of a cycle, so 
we remove it and begin searching again with Ok-i (if k = 0, we choose a new oq 
arbitrarily from among the remaining vertices of V). We continue this process until 
we have removed all the vertices of V or we have found a cycle. 

In the former case, if we list the elements of V in the order in which they were 
removed from V, we have some ordering oq, ai, a2, . . . , a^ in which i < j implies 
ai 7^ Oj (since Oi could otherwise not have been removed before Oj). Thus, we can 
define an acyclic orientation A on Gr by simply saying that whenever {ai,aj) is 
an edge in Gr, Oj -<a Oi. Let a be the tiling for which GR{a) = A. Then if ic,i is 
the last tile that occurs in our ordering, we know that tc^iia) contains the leftmost 
squares of rows c, c-t-1,..., c + n— 1. We know that the next to last tile td,j{a) 
contains the leftmost remaining squares of rows d,d+l, . . . ,d + n — l, and so forth. 
Thus, given our ordering, it is easy to explicitly determine all the tiles of a in linear 
time. Note also that this algorithm yields an explicit proof of the fact (mentioned 
earlier) that every acyclic partial ordering of a graph has an acyclic extension. 



4. General Tilings And Abelian Height-Functions 

4.1. Examples. The constructions in this paper allow us in principle to work out 
the maximally informative abelian height function for any polyomino tiling problem. 
The following facts, for example, are simple to verify: 

1. Given any coloring of the squares of R, a generalized ribbon tiling is a tiling 
in which each tile contains exactly one square of each color. Suppose, for 
example, we color the square at position (x, y) one of mn colors according to 
the values of x (mod m) and y (mod n) . In this case, the canonical abelian 
height space is Z™""-"^. (When m — n = 2, this is the square/skew-tetromino 
tiling problem mentioned in |l7|.) 

2. Let T be the set of all polyominos with two black and two white squares in the 
usual chess board coloring. (These are squares, skew-tetrominos, and length- 
four vertical and horizontal bars.) The canonical abelian height function space 
is Z X Z/2Z. 

3. If T contains all trominos, then the canonical abelian height space is trivial. 

4. Let T be the set of all rectangular m x 1 tiles and 1 x n tiles. Then the 
canonical abelian height space is Z^™"^^*^""^^. (This problem was studied in 
P], where connectedness under local moves was proved using a non-abelian 
height function.) 

5. If T contains only horizontal dominos, then the canonical abelian height space 
is infinite dimensional. 
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Figure 7. Three possible spanning trees for a single tile. 

4.2. Defining height functions at all vertices. In general, the height function 
ha is only defined at the vertices of R that are on the boundary of tiles of a. One 
(admittedly awkward) way around this is to replace each tile with a spanning tree 
of the squares in that tile, and consider an edge in t to be on the boundary if the 
squares it divides are not adjacent as members of the tree. (See Figure 0.) That is, 
if vi and V2 are adjacent vertices in i?, we say ha(v2) — ha{vi) is equal to (f){vi,V2) 
when the squares incident to (i'i,W2) are either contained in different tiles or are 
contained in the same tile t but are not adjacent as members of the spanning tree 
of t; this is enough to determine the value of ha at every vertex v in R — and this 
is now true even if we allow T to contain non-simply-connected tiles. 

Of course, the height values now depend on the particular choice of spanning 
tree for each tile; one might even consider each way of choosing a spanning tree of 
a tile i to be a different tile in its own right. In this perspective, T becomes a set 
of tile-spanning-tree pairs. Changing the spanning tree corresponding to a given 
tile is a sort of local replacement move. 

4.3. Can distinct tilings have the same height function? We saw in the 

examples that the abelian height function space for tromino-tilings is trivial; in this 
case, of course, all tilings of R by trominos have the same (trivial) height function 
mapping all vertices to zero. 

We will now describe a condition on tile sets T that will ensure that two distinct 
tilings do not have the same height function. Given a map (/) from the squares of R 
to the abelian height group G, we say a tile t is admissible by (f) if (/)(t) = 0. We say 
t is irreducible with respect to (j) if there is no way to partition t into two or more 
smaller polyominos each of which is admissible. A tile set T is irreducible if each 
tile is irreducible with respect to the map (f> from the squares of R to the canonical 
abelian height space determined by those tiles. Note that these definitions depend 
only on the (canonically determined) values of 4> on squares s; they do not depend 
on any arbitrary choices we made when choosing the value of on individual edges. 

Now, consider two tilings a and P oi R with the same height function h. We 
claim that each connected component u of the intersection of any tile in a with a 
tile in /3 satisfies 0(m) = 0. To see this, let vo,vi, . . . , Vk = vq be the vertices of a 
cycle surrounding u. Since each edge in this cycle is a boundary edge of either a 
or /3, it follows that 



^(u) 



E' 



y^ 4){vi,Vi-i) = ^ h{vi) - h{vi^i) = 0, 



z=l 



z=l 
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where the second sum is over clockwise edges on the boundary of u. Unless u is 
a tile in both (3 and a, it must be contained in a larger tile t in one of the two 
tilings, say (3. All the connected components of intersections of t with tiles of a are 
admissible, and it follows that t is reducible with respect to 0. We can now deduce 
the following: 

Proposition 4.1. If T is irreducible, then no two tilings of R by tiles in T have 
the same height function. 

This does not mean that height functions for reducible tile sets cannot provide 
useful information, but if a tile set is reducible, there may be aspects of the tiling 
that arc not captured by the height function. A good example is the tile set T 
containing all domino tiles and all skew tetromino tiles. The abclian height function 
for tilings by these tiles is simply the usual domino-tiling height function, where a 
skew tetromino is treated as a pair of dominos. 

4.4. Which functions are height functions of a tiling? Assume our height 
functions are defined at every vertex. Then for each oriented edge e — (i;i,W2) we 
define Ag to be the set of all differences ha{v2) — ha(vi) that can possibly arise for 
a tiling a of any region R. Recall that: 

where the latter sum ranges over all squares passed on the left by the edge (wi, U2). 
Thus the set Ag is equal to (t>{vi, V2) + Lg where Lg is the set containing and all 
possible square sums ^ (f>{s) over sets of squares that can occur in the left side of 
a tile divided by e. Clearly, Ag == — Ag-i. 

We say a function h from the vertices of i? to G has proper differences with 
respect to the tile set if 

1. h{v2) — h{vi) = (f>{vi,V2) whenever (wi,f2) is a boundary edge of R. 

2. h{v2) — h{vi) e A(^,j^ „2) whenever (vi, V2) is an interior edge of R. 

3. h{vo) = for some reference vertex vq on the boundary of R. 

We say a tile set T is complete if it is irreducible and for every simply connected 
region R, every function on R with proper differences is the height function of 
a tiling of R by tiles in T. For this definition, cj) is assumed to be the usual 
map from edges into a group G containing the canonical abelian height space as a 
subgroup. Like the definition of irreducibility of T, this definition is independent 
of any arbitrary choices we made in choosing (j>. To see this, note that if we add 
to a function (j)'{e) on the edges of Z^ such that ^'(s) = for any square, this 
simply amounts to adding a tiling-independent function to each ha and adjusting 
the Ag values accordingly. When a tile set is complete, there is a precise one-to-one 
correspondence between tilings of R by that tile set and height functions on the 
vertices of R with proper differences. 

The fact that dominos are a complete set of tiles has been very important to the 
study of domino tilings; it is often easier to work with the space of height functions 
satisfying difference restrictions than to work directly with domino tilings. Com- 
pleteness of domino tilings implies that existence algorithms and random sampling 
algorithms for domino tilings can all be derived from corresponding algorithms for 
spaces of height functions with proper differences. 
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4.5. Generalized Ribbon Tiling Completeness. We show here that any set of 
generalized ribbon tilings is complete. Given a coloring of the squares of T? with n 
colors, recall that a generalized ribbon tile set T consists of all polyominos (together 
with corresponding spanning trees) that contain exactly one square of each color. 

Let G be the canonical abelian height space for the generalized ribbon tile set 
T. Although we do not in general know the structure of G, we can determine the 
structure of iJ, the quotient of G with the relations (/'(si) — (i>{si) — whenever s\ 
and Si have the same color. Let ^ be the composition of ^ with the quotient map 
sending G to H and observe that the value of ^p{s) depends only on the color of s 
(the same is not necessarily true of 0(s)). Let Ci be the value of ip{s) when s has 
color i. 

Now H is precisely the group generated by the e^ subject to the one relation 
SfcTo ^i — 0- (Because this relation is enough to force tp{t) to be zero for each 
t G T, no additional relations are required.) It is easy to see that these tiles are 
irreducible: no nonempty proper subset of the e^ sums to zero in H, so each sum 
of (j){s) over a nonempty proper subset of the squares of a tile has nonzero image 
in H and is hence a nonzero element of G. 

Now, to show that T is complete, we must show that every h mapping the 
vertices of i? to G with proper differences is the height function of a tiling a. 
Suppose h has proper differences. Given an edge e = (wi,i'2) in R, we define 
6(e) = h{v2) — h(v\) — (/'(e). By the definition of proper differences, 6(e) is the sum 
of 0(s) over some set S'(e) of squares s — and there exists at least one tile t E T 
such that e crosses t and passes those squares on the left. In particular, this implies 
that ^(e) contains at most one square of each color and cannot contain a square of 
every color. Put differently, the set G(e) of colors that occur in squares of S{e) is a 



iec(e) ^3 



proper subset of {1,2,... , n}, and the image 6(e) of 6(e) in H is simply J2 

We say e has color j on its left if j € G(e). Note that if e has a square of color j 

immediately to its right and h has proper differences, then G(e) can never contain 

J- 

Now we will construct the tiling a for which h is the height function. Let ai, 

a2, 03, and 04 be the clockwise oriented edges of a square s with color i. If /i is 
a function with proper differences and wo,i'i,^'2ji^3, ^4 — vq are a path of vertices 
surrounding s clockwise, we have: 

44 4 

= ^ h{vk) ~ h{vk-i) = ^ (j){ak) + b{ak) = 4){s) + ^ b{ak) 

k=l fe=l fe=l 

Taking the image of the latter sum in H gives 

4 
0^e, + J2Kak) = ei + Y^ Y^ e^ 

fc=ljeC(afc) 

Since none of the C{ak) can contain i; it is easy to see that the only way ^ 6(0^) 
can be — e^ in _ff is if the sum includes exactly one term of type ej for each j ^ i- 
Thus, for 1 < k < A, the G{ak) are disjoint, and their union contains all colors 
except i. 

Now, we say si and S2 arc tile adjacent if they are adjacent and the edge e between 
them has a nonempty S{e). We define our tiles to be the connected components 
of the tile-adjacency relation. If we can show that every tile contains exactly one 
square of each color, we will have proved the completeness of T. 
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Suppose So has color ?, and choose some j ^ i. We aim to show that sq is in 
the same tile as a square of color j. Observe that at least one edge of the four 
clockwise edges of S has color j on its left. Let si be the square adjacent to Sq that 
is incident to the edge e for which this is the case. Clearly, si does not have color i. 
If it has color j, we are done. Otherwise, suppose ai, 02, and 03 are the three other 
clockwise oriented edges of si. We know that the union of the C{ai) and C(e~^) 
contains all colors except the color of si. Since C(e~^) U C(e) = {1,2, .. . ,n}, it 
follows that each of the C{ai) is a proper subset of C(e) and that one of them 
contains j. Now, we let S2 be the square adjacent to Si and incident to the edge 
for which this is the case and repeat this process. Since the proper subsets cannot 
decrease in size indefinitely, we must eventually reach a square of color j. A similar 
argument shows that such a path starting at sg never reaches another square of 
color i. It follows that each tile contains exactly one square of each color and that 
the adjacency relation describes a spanning tree of that tile. 

Using similar arguments, one can show that the tile set containing all 1 x fc and 
m X 1 rectangles is complete. 

4.6. Completion of a Tile Set. A good example of a tile set that is irreducible 
but not complete is, when n > 3, the set of ribbon tiles of order n with all but one 
of the possible ribbon tile shapes. 

When is it possible to add tiles to an incomplete tile set T to make it complete? 
First, it is clear that we will have to add as a tile every finite region R on which 
there exists a height function h with proper differences such that the set S of edges 
(w 1 , W2 ) in i? for which h(v2) — h{v2) = 4>{vi,V2) does not divide the squares of R into 
more than one component. (There may be an infinite number of such regions, and 
they may be arbitrarily large.) Note that since all of these new tiles are already 
admissible by 0, adding these tiles will not change the canonical abelian height 
space, and we can use the same 4> to define height functions for the new tiling. 

If each of these new tiles is irreducible, the result is a complete set of tiles, called 
the completion of T. 

4.7. Non-Existence Proofs for Complete Tile Sets. Conway-Lagarias M and 
Thurston p7[ | observed that a simple necessary condition for the existence of a tiling 
on R is that the height function be well-defined on the boundary of R. That is, the 
group-wise sum of (/)(e) around the all the edges e oi R must be the identity. 

It is not difficult to describe much stronger necessary conditions. If a tile set T is 
complete and i? is a simply-connected region, then (assuming the previous condition 
is satisfied) the question of whether R can be tiled by T reduces to the question 
of whether there is a height function on the vertices of R with proper differences. 
Even if T is not complete, the existence such a function is a necessary condition 
for the existence of a tiling of R. Thus, the following interpolation problem is of 
interest: 

Given a simply connected region R, a function h defined on all boundary vertices 
of R, and a set Ag for each oriented edge (with Ag = — Ag-i for all e), is it possible 
to interpolate h to all vertices of R in such a way that h{v2) — h{vi) G A(-„j „,) 
whenever (t^i, V2) is an edge in Rl 

When the height function space is a rank n free-abelian group, this is a special 
case of integer programming, which is in general NP-complete. But suppose we 
embed the height space in R" and let Ae be the convex hull of the set of points 
in Ae. Then a relaxed version of the problem asks for an interpolation of h to 
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the interior squares of R where h is aUowed to assume any values in R" and the 
differences are merely required to lie in Ag. Since each Ae is a convex polyhedron, 
this is now a linear programming problem that can be solved in polynomial time. 
If we prove that there is no solution to the relaxed interpolation problem, we know 
R cannot be tiled by T. 

4.8. Note on Multiply Connected R. Suppose that R is simply connected ex- 
cept for k "islands," labeled Ii, . . .Ik] we assume that the boundary edge/ vertex 
sets of the boundary of each Ik and the outer boundary of R are disjoint, but that 
each li is connected. 

The problem with defining height functions on R is that as we traverse the 
boundary of a given Ik clockwise and come back to the same place, the height 
function changes by (f>{Ii); thus, if (f>{Ii) 7^ 0, then the height function ha is not 
well-defined as a single-valued function. There are two ways around this problem. 
One is simply to accept that ha is a multiple-valued function; only the differential 
\/ha{vi, V2) = ha(v2) — ha{vi), is well defined, and this can be integrated to produce 
a single- valued ha on the vertices of any simply connected subset of R — but not all 
of i?. 

Second, we could choose some replacement (f>' for (f> satisfying: 

1. (/)'(/i) = for each island /j. 

2. (/)'(s) — <j){s) for each s in R. 

It is not hard to see that such a (j)' always exists. If we define ha using (j)' instead 
of (j), then it is everywhere defined. We require ha{vo) = for some outer boundary 
vertex vq, and we also choose one vertex Vi on each island i. Clearly, the value of 
ha on all boundary vertices is determined by its value on the Vi. 

In the special case of ribbon tilings, we define a simple local replacement move 
on _R to be a move that involves removing two tiles ii and ^2 that together comprise 
a simply connected two-tile region R' and replacing them with another pair of tiles. 
When R is simply connected, no tiling of R ever contains a pair ii and ^2 of tiles 
whose union is connected but not simply connected — otherwise, there would be a 
hole between them that could not be filled by order-n ribbon tiles. Thus, when R 
is simply connected, every possible local replacement move is simple. It is not hard 
to show in general that two order-n ribbon tilings a and (3 on R can be connected 
by a sequence of simple local replacement moves if and only if ha{vi) = hi3{vi) for 
each i. The reader may also check that when n < 4, every local replacement move 
is simple. 

5. Open Problems 

The domino tiling problem is perhaps the most well understood of all tiling 
problems. Height functions and local replacement moves have played a fundamental 
role in the proofs of many domino tiling theorems, including the following. 

1. Proofs of local- move connectedness. 

2. Algorithms for computation of minimal- length local move sequences connect- 
ing two tilings. 

3. Linear algorithms for determining when a region R has a tiling. 

4. Mixing time bounds on random processes based on local moves. |12], |gg 

5. Perfect sampling algorithms based on "coupling from the past." |1S| 
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6. Asymptotic entropy computations for random tilings and descriptions of finite 
and asymptotic Gibbs measures. |^, g 

7. Law of large numbers and large deviations principles for asymptotic height 
function shapes of random tilings. S 

8. Conformal invariance of certain tiling properties in the asymptotic scaling 
hmit. |ll| 

These results and references are ridiculously far from exhaustive; there have been 
numerous papers relating to each of the above topics and many others. We include 
the references not to assign credit, necessarily, but to provide pointers to papers 
with more complete bibliographies. Because ribbon tilings are in many ways the 
simplest polyomino tilings that generalize domino tilings, the following question is 
very natural: 

Question 5.1. Which of the above results can be generalized to ribbon tilings? 

An additional task is to generalize to other tiling spaces with abelian height 
functions. 

Question 5.2. Can connectedness under local replacement moves be generalized to 
the other tiling problems in this paper? In particular, is there a local move result 
for every complete, translation-invariant tile set? 

Finally, the author hopes this paper will be a step towards a more thorough 
understanding of at least the abelian case of the Conway-Lagarias height function 
theory. The following is a more theoretical question: 

Question 5.3. Is there a simple classification of the complete, translation invari- 
ant tile sets with finitely many tile shapes? Given such a tile set T, is there always 
a polynomial algorithm for determining whether a simply connected region can be 
tiled by T ? 
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